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Abstract. We perform the analysis of the trispectrum of curvature perturbations 
generated by the interactions characterizing a general theory of single-field inflation 
obtained by effective field theory methods. We find that curvature-generated 
interaction terms, which can in general give an important contribution to the amplitude 
of the four-point function, show some new distinctive features in the form of their 
trispectrum shape-function. These interesting interactions are invariant under some 
recently proposed symmetries of the general theory and, as shown explicitly, do allow 
for a large value of the trispectrum. 
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1. Introduction 

Inflation [1] is one of the central pillars of modern cosmology. Not only it provides 
a natural solution to the flatness, horizon and monopole problems of standard Big- 
Bang cosmology, but can also explain the production of density perturbations in the 
early Universe which then lead to LSS [2, 3, 4, 5, 6] in the distribution of galaxies and 
temperature anisotropies in the CMB [7, 8, 9, 10, 11, 12]. 

Besides the simplest single-field slow roll inflationary model, many other inflationary 
mechanisms have been proposed since inflation was first introduced, and all are 
compatible with the CMB and LSS observations. In order to probe deeper into the 
dynamics of inflation and to remove the degeneracies generated by the many models so 
far proposed, one might study observable quantities which are sensitive to deviations 
from Gaussianity [13]: starting from the three-point function [14, 15], one then considers 
the trispectrum [16] and in general higher-order correlators as well as loop effects 
in the power spectrum [17]. Such non-Gaussian features will depend on the various 
interactions characterizing any given inflationary model in the form of self-interactions 
of the inflaton, its coupling with gravity and interactions with other fields in the case of 
multi-field inflation (see, e.g., [18, 19, 20] for comprehensive and updated reviews). These 
investigations are spurred by the fact that the continued analysis of WMAP data [12] 
and the recent launch of the Planck satellite [21, 22] provide the exciting opportunity 
to actually test the predictions of this zoo of models at the level of bispectrum and 
trispectrum of curvature perturbations. A very useful tool in analyzing the possible 
signatures of the different inflationary models is given by the effective field theory 
approach to inflation recently introduced in [23] and further expanded in [24, 25, 26, 27]. 
There are various advantages in employing this formalism. Indeed, it provides a unifying 
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perspective on inflation in that it automatically accounts for many known inflationary 
mechanisms. To each set of interactions for a given inflationary model there corresponds 
in the effective Lagrangian a linear combination of operators obtained by turning on 
and off some specific coefficients that regulate the weight of the operators (we introduce 
these coeffcients later on and refer to them as M„'s). The unifying power of the effective 
field theory approach is quite manifest in that, in principle, it allows these coefficients 
considerable more freedom than what they are granted in any specific inflationary model. 
In fact, by being for the most part free parameters (a couple of these coefficients are 
to obey some inequalities if one wants, as we do, the generalized speed of sound to 
be smaller than unity), the M n 's allow for the description of known interactions with 
relative weights which would otherwise be fixed, so by employing effective field theory 
one enlarges the region of the parameters space than can be spanned. Besides that, in 
the effective Lagrangian some of the M n 's multiply curvature-generated operators (see 
below Eq. (2)) that are sometimes neglected but should in principle be studied as, in 
fact, their contribution can be relevant [28] and increase the dimension of the parameters 
space of the theory. 

In this paper we use effective field theory techniques to study the trispectrum of a 
very general theory for single field inflation. In particular, we concentrate on the 
contributions of many novel curvature-generated terms as well as interactions that 
characterize Ghost inflation [29]. As an ordering principle among the numerous 
interactions such a general approach comprises, we employ two additional symmetries 
of the action recently introduced in [30, 31] (see also [32]) and examine the cases when 
one or both of these requirements are imposed on the theory. Following [33] , we analyze 
the shapes of these terms in four different configurations so as to identify distinctive 
effects in the trispectrum shape-function from the various interactions. We first analyze 
the contribution to the scalar exchange diagram due to a curvature-related term which 
generates an interesting flat shape for the bispectrum [28] . This term produces a shape 
function which is, in some configurations, different from all the shapes due to leading 
interactions in general single-field inflation models [33]. f We calculate and plot the 
contributions to the contact interaction diagram by several interactions: we rediscover 
the shape-function due to the leading fourth-order interaction in ghost inflation first 
obtained in [30] and plot it in new configurations; an analysis of other interesting term 
is also performed allowing us to identify novel distinctive features of curvature-related 
interactions. 

The paper is organized as follows. In section 2 we build on [24] to introduce the very 
general effective theory we employ in all subsequent calculations. The reader who is 
familiar with this procedure might want to skip this part and start directly from Eq. (9). 
In section 3 we characterize the various interaction terms of the theory according to 
their behaviour under the action of two specific symmetries. Afterwards, in section 4, 
we proceed by briefly outlining the tools of the IN-IN formalism successively employed 

| In these models the inflaton Lagrangian is an arbitrary function of the inflaton and its first derivative. 
The theories we will consider here further generalize these models. 
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in the trispectrum calculations. A separate analysis of the scalar exchange and contact 
interaction diagrams contributions to the trispectrum is performed. Section 5 contains 
a summary of the findings and comments on further work. In Appendix A we report 
the details of the scalar exchange diagram calculations. In Appendix B we show the 
reliability of our simplifying assumption on the classical solution to the equations of 
motion for the general theory considered in this paper. 

For the sake of clarity we stress here that whenever we refer to general single- 
field inflation models, often when elaborating on the results of [33], we are dealing 
with theories which account for a great fraction of known inflationary models (DBI, 
K-inflation etc.), but still miss an important subset including Ghost inflation and 
curvature-generated interactions in general. On the other hand, the effective approach 
of [24] we employ here also covers the latter. 



2. The Hamiltonian up to fourth order 

We will follow the effective theory approach first introduced in Ref. [24] in order to 
write down the complete theory of single-field models of inflation up to fourth order 
in perturbations, subjected to the sole requirement of an approximate shift symmetry, 
it — > 7i + const, on the scalar degree of freedom tt. Let us give a brief account on how 
to obtain the main formulas. We start with the scalar field responsible for inflation, 
which is split as an unperturbed part plus the fluctuation: 

<j)(x,t) = Mt) + s<j)(s,t). (i) 

For reasons that will soon become clear, one chooses here to work in the the comoving 
(or unitary) gauge for which 5<p = [23]. As a result, the action will no longer 
be invariant under full spacetime diffeomorphisms (diffs) but only under the spatial 
reparametrizations. This is the starting point to write the most general unitary gauge 
space diffs invariant Lagrangian at the desired order in perturbation theory [24]: 

S = J d A x^F(R, upa ,g 00 ,K, u ,V„t), (2) 

where is the extrinsic curvature and the indices on the metric entries g 00 are free 
indices. Taking into account the fluctuations around a FRW background, one obtains 
the following action 

S= J d 4 x^/=g~ -M^R + M^Hg 00 - (3H 2 + Hj + 

J2 F (n \g 00 + 1, 8K^ 5R, upa , V M ; t)] , (3) 

n>2 

where the fluctuations contained in are at least second order. The next step is 
restoring full spacetime diffs invariance. To see how it works, we take from [24] the 
following sample terms in the action: 
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d 4 Xv ^g-[A(t) + B(t)g 00 (x)} . (4) 
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Consider the time reparametrization: t — >■ t + £0(2^); % — > %\ under its action (and 
after a simple variable redefinition) Eq. (4) reads: 



(5) 



At this stage the procedure we will adopt consists in promoting £ to a field, $,q{x) = 
— tt(x) and requiring the following gauge transformation rule: tt(x) — >■ tt(x) — £,o( x ) on 
71. With this assumption in place, the above action is invariant under full spacetime 
diffeomorphisms. The scalar degree of freedom it makes its appearance in the time 
dependence of the A(t + tt), B(t + it) coefficients and in the transformed metric. This 
procedure is essentially the same as the one of standard gauge theory: a Goldstone 
boson which transforms non-linearly under the gauge transformation provides the 
longitudinal component of a massive gauge boson. For high enough energies, the 
Goldstone becomes the only relevant degree of freedom. This is the so-called equivalence 
theorem. The same is true for our case: for sufficiently high energy the mixing with 
gravity becomes irrelevant and the scalar tt becomes the only relevant mode in the 
dynamics (decoupling regime). One then needs to identify the scale of the energy above 
which this approximation holds (on the other side of the energy range one always keeps 
in mind the upper energy threshold that comes with the use of effective theory). This 
procedure puts some bounds on the values of some of the coefficients that drive quadratic 
operators in the action [24, 28]. We will not be concerned with these issues because most 
of the interaction we will be analyzing start at third order in perturbations. It suffices 
here to say that, since one is concerned with correlators just after horizon crossing, the 
decoupling procedure works as long as the decoupling energy is smaller than the Hubble 
rate H . 

From now on we will work in the decoupling regime. In considering the terms of Eq. (3), 
we will therefore use only the unperturbed entries of the metric tensor. To write the 
effective Lagrangian up to third order, we start from Eq. (3) and follow the algorithm 
given in [24]. Fluctuations are encoded in the terms. In order to be as general as 
possible, one must also include all possible contributions coming from extrinsic curvature 
K^ v terms. 

Following the procedured outlined above, the third and fourth-order Lagrangian is 
obtained. We will use the third-order expression to calculate the contribution to the 
trispectrum of curvature perturbations arising from the scalar exchange diagram [28]: 
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The above action and its extension at higher perturbative orders covers many 
inflationary theories, providing a unifying perspective which is hard to obtain without 
an effective approach. Indeed, by switching on and off a single M n = M n , M n coefficient 
one has control over all corresponding operators in the action (the bar on the M n 's 
signals that these are curvature-generated terms). The hope is to be able to identify 
distinctive features for as many as possible different combinations of the M„'s in the form 
of specific patterns they produce in the shapes of the various correlators of curvature 
perturbations. The degeneracies among the results for different inflationary mechanisms 
that will inevitably arise might be removed by a joint analysis of the different n-point 
functions, starting with the bispectrum, the trispectrum, loop corrections to the power 
spectrum and so on. Let us briefly go through some of the main features of the third 
order effective action above. All the comments can be straightforwardly extended to 
the fourth-order expression as well. 

• Consider only the quadratic terms: for M 2 = M 1)2 ,3 = one recovers the usual 
quadratic Lagrangian for the fluctuations, with sound speed (? s — 1. Switching 
on M 2 corresponds to allowing models with sound speed smaller than unity, 
\j(? s = 1 — 2M$/(Mp l H), which are often linked to a high level of primordial 
non-Gaussianity [35, 24] as for DBI inflation. Further allowing for a non-zero M 2j3 
in the de Sitter limit, one recovers Ghost Inflation [29]. Similarly, having all M's 
set to zero, and going to third and higher order with the M's, one retrieves the 
interactions characterizing DBI inflation [34, 35, 33] and K-inflation [36, 37] theories 
and others. § 

• The action in Eq. (6) contains in principle additional terms but, being interested in 
those generating large non-Gaussianities, a selection has been made. Specifically, 
at every order in fluctuations and for each M n coefficient, only leading terms are 
considered. One starts from the realization that, even for the most generic quadratic 
action the following estimates hold at horizon crossing [28]: 



where c s is a sort of generalized speed of sound. 

Also, as mentioned before, in the action the scalar appears only trough its 
derivatives. When faced with a given M n multiplying terms of the same 
perturbative order with a given number of derivatives one therefore knows that 
the leading term will be the one with the most spatial derivatives. 
There is also the comparison between the same perturbative order but different 

§ See [38, 39] for some recently introduced examples that require, in order to be recovered in the 
effective field theory approach, that one relaxes the implicit assumption of a shift symmetry for the 
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M n terms to be made. All non-zero coefficients in front of the various operators 
might be assumed to be of the same order [29]; interestingly, in [27], employing 
renormalization and unitarity arguments, a natural (relative) value of 1/c 2 . was 
obtained for the M|, M| coefficients. In general, we shall not restrict ourselves to 
these situations. Considering theories with a speed of sound different from unity 
and allowing for M n 's of different orders greatly increases the number of viable 
terms for large non-Gaussianities. Let us consider an illustrative example. Take 
the interaction terms 

WW 4,M 2 . M 7 (t) (dfnf 

At horizon crossing, the region from which we expect the main contribution to 
n-point functions, the comparison reads like M^HH 2 /^ 2 -H- M 7 H 6 /c^. In a 
Lorentz invariant theory with coefficients of the same order the first term would 
clearly prevail. Allowing c s <C 1 makes the comparison less obvious and an 
Mi 3> M 2 further strengthens this point. Here a word of caution is in order: 
from simple dimensional analysis a term with ever increasing spatial derivatives 
will have an M n coefficient with smaller and smaller exponent (counterbalanced 
by an higher exponent for H) and must eventually be subleading with respect to 
the contributions with fewer derivatives. This is because in the effective theory 
approach one is roughly making an H/M expansion (with M being the scale of the 
underlying theory) and, although one can fully employ the freedom to have M n 's 
of different size up to some perturbative order in order to resuscitate interesting 
contributions to the correlators, the M n -driven contribution must eventually (from 
some n on) cease to be relevant. 

Employing the same calculational algorythm first introduced in [24] and used in writing 
the complete third order action above, we obtain the most general fourth-order action 
in this set up: 



<S 4 = / d 4 x\/^g 



^,(,)^ + 2M , ( ,)«W + ^( t )V 
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2 V 2a 6 a 6 
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Note that, as pointed out in [40], starting at fourth order in perturbations, one 
cannot immediately read off the Hamiltonian from the expression of the Lagrangian, in 
other words H = —L does not hold here. We use the results one obtains by adopting 
the correct procedure which was outlined in detail in [40]. 

Let us split the interaction Hamiltonian we will be concerned with as H int = H 3 + i7 4 ; 
one can prove that the overall interaction Hamiltonian is then: 

-L3-L4 + J Sx^f- f 1 



H, 



int 



2M 2 4 + H 2 eM P - 3HMf 
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(9) 



where the above terms besides — (L 3 + L 4 ) are all at fourth order in perturbations. 



3. Symmetries 

Having written the complete Hamiltonian, we now proceed to calculate the four-point 
function contributions arising from interaction terms at third and fourth order. We 
employ here the IN-IN formalism [42, 43, 44, 45] and conveniently split the contributions 
to the four-point function as the ones arising from terms that make up the contact 



8 



interaction diagram and the ones that generate the scalar exchange diagram as in the 
figure below. 

X 

Figure A: On the left, the scalar exchange diagram. Contact interaction diagram on the right 

It is useful at this stage to offer some comments on the calculations we are going 
to present. As mentioned, the literature already contains a thorough analysis of the 
trispectra for general single-field inflation models, see for example [33]. Work on the 
four-point function for ghost inflationary models has recently been presented [32, 30]. 
Our starting point, being based on a comprehensive effective theory, clearly encompasses 
all these models. Working with the effective Hamiltonian above translates into many 
immediate advantages as listed before but, on the other hand, in calculating the resulting 
four point function, one faces a substantial number of terms and it is therefore natural 
to look for some ordering principle which would single out some contributions to the 
trispectrum as the leading ones and allow us to concentrate on them only. In this context 
employing a symmetry for the whole theory can prove very useful. Indeed in [31, 30] 
the authors consider only those allowed by a particular (approximate in [31]) symmetry 
of the action, respectively: 

SI : tt —> —it ; S2 : vr ->• -vr and t -t . (10) 

We plan here to employ our general effective theory to show that, allowing some freedom 
on the M n coefficients that modulate the various terms in the third and fourth order 
action, within each one of the two distinct and quite restrictive symmetry requirements 
above there are novel curvature-generated terms in the action that should not be disre- 
garded as negligible and that, furthermore, show some distinctive features in the shapes 
of the trispectrum. We will also describe terms allowed by both the symmetries in 
Eq. (10) combined. Of course, one need not employ symmetries to switch on or off any 
specific operator in the action. Most of the contributions are indeed freely adjustable 
by the correspondent M n coefficient, a procedure which is, in principle, legitimate since 
the underlying theory is unknown. We choose here to restrict ourselves to considering 
only symmetry-abiding terms. Let us comment on each one of the symmetries. 
SI is built upon the following considerations. Often the same M n coefficients multiply 
terms of different perturbative orders; consequently the amplitude f^L of the 3-point 
function will be related to the amplitude of higher order correlators, notably to t^l, the 
amplitude for the four-point function. Whenever the leading part of the trispectrum is 
generated by these types of M n 's one can estimate that for its effect to be observable 
t nl has to be five orders of magnitude larger than f NL [31], which leaves little room for 
feasible models. On the other hand, one quickly realizes those M n 's whose first term 
starts only at the fourth perturbative order (M 4 , Mi ...in Eq. (8)) are not plagued by this 
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problem. This then represents a natural way to obtain inflationary models which allow 
a large, detectable trispectrum untied to the interactions which make up the bispectrum 
(which might well be small now).|| Indeed, in [31] the authors investigate on the size 
of all the interactions driven by the M 4 , Mn, Ni % coefficients in Eq. (8) and show that 
the leading interactions driven by these parameters are all consistent with the tt — > —n 
prescription and are expected to give a comparable signal + . By construction then, the 
terms in the interaction Hamiltonian that are going to contribute to the trispectrum 
and be consistent with the reasoning that inspired the SI symmetry are only some of 
the ones that will make up the contact interaction diagram, namely those whose lowest 
order interaction is already at fourth order. This limits us to the contributions regulated 
by the following coefficients: M 4 , M W ..M 15 , Aq..7V 5 . 

S2 symmetry, on the other hand, does not prohibit third order interactions, indeed in 
[30] the interaction 7r(V7r) 2 is considered and, by inspection of Eq. (6), one can see that 
also other terms are allowed, the one regulated by M 5 and, notably, the M & ii(dijii) 2 / a A 
term. The M 6 -driven term is particularly interesting because its contribution to the 
bispectrum calculations of [28] generates an interesting flat shape. The scalar exchange 
diagram will then be built out of the third order S2-obeying terms in the action. In 
particular, inspired by previous findings, we are going to give a detailed account of the 
Mq contribution. 

If both SI and S2 are to be enforced one must also exclude from the list of Sl-abiding 
interactions the ones multiplied by M±o, M 13 , M 14 , M 15 . A more clear picture of the 
situation concerning the various symmetries is presented in Table 1 below. 



Table 1 



Coefficients 


M 2 


M 3 


M 4 


Mi 


M 2 


M 3 


M 4 


M 5 


M 6 


M 7 


M 8 


M 9 


SI 


X 


X 


/ 


X 


X 


X 


X 


X 


X 


X 


X 


X 


S2 


/ 


/ 


/ 


X 


X 


X 


X 


/ 


/ 


X 


X 


X 


Coefficients 


Mio 


Mn 


M 12 


M 13 


Mia 


M 15 


Ni 


N 2 


N 3 


N A 


N 5 


/ 


SI 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 




S2 


X 


/ 


/ 


X 


X 


X 


/ 


/ 


/ 


/ 


/ 





The Coefficients marked with " /" in correspondence of a given symmetry S are 
S-invariant, those marked with "X" violate the S symmetry. 



|| One needs also to check that the interactions driven by coefficients that multiply also third order 
fluctuations do not become important in the form radiative corrections to the bispectrum. This check 
is done in [31] and ensures that loop corrections of those terms are not relevant. 

In the same spirit of the analysis done in [28] for all curvature-generated terms at third order, the 
authors of [31] consider in the v2 of their paper some extrinsic-curvature terms generated at fourth 
order. They also comment on their importance in near de Sitter limit and their conclusions apply to 
our Mu,Ni parameters. 

+ It would be interesting to understand to what kind of models, in terms of the fundamental scalar 
field, the simple resulting effective Lagrangian corresponds in this case. 
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Note that each M n coefficient might multiply many interactions at each 
perturbative orders and therefore we mark the coefficient as invariant under a symmetry 
when all the leading interactions it multiplies are invariant under SI or S2. Determing 
the properties of the coefficients in the second row requires no effort, as one can easily 
verify these M n 's first appear in the action as multipliers of fourth-order terms. Things 
are less linear with the coefficients in the first row (except for M 4 ) as they appear at 
fourth order both multiplying bare interaction terms and multiplying other coefficients 
as well as interaction terms (for an example of the latter case see the terms written 
explicitly in Eq. (9)). They also appear at third and some also at second order in 
perturbations. One then must carefully check that, given a particular coefficient M n , in 
none of the interactions it multiplies at any order the leading terms violate the symmetry. 
For M2, M3, M5, Mq in the first row one can verify after some checks that these terms all 
parametrize indeed approximately invariant interactions upon requiring the coefficient 
M| to be much smaller than the typical M n such as M 2 ..M 6 . This is because in the 
fourth-order Hamiltonian in Eq. (9) there are terms of the form 



which one then assumes to be subleading. We stress this point because it emerges clearly 
and naturally in the effective theory approach. 

4. Trispectrum 

4-1. IN-IN Formalism 

We are going to employ the IN-IN formalism to calculate the four point function of 
curvature perturbation. The most general and compact expression for such a quantity 



(nic*iCwC«c*4(f)in> = (o|T{e i /^ d3 ^'«W}a 1 a 2 a3a4(t)T{ e -^ d3 ^ d '" w W}|o) 

(12) 



where T and T indicate respectively anti-time order and time order operations, |0) and 
\fl) stand for the vacuum of the free and interacting theory. 

Expanding both the exponentials in Eq. (12), we single out the first non vanishing terms 
that will contribute to the scalar exchange and contact interaction diagrams. 




x {M 2 4 , M 3 4 , Mf , M 6 2 } x (S2 - violating interaction) 



(11) 



is: 



(fl|CfciCfc2C*3G4(t)|n> = 



(0\f{i / d 3 xdtH 3 (x)}( kl ( k2 ( k3 ( k4 (t)T{-i / d 3 x'dt"n 3 (x')}\0) 




+ (0\f{ 1 - [ I d 3 xdt d 3 x dt'H 3 (x)H 3 (x')}( kl ( k2 ( k3 ( M (t)\0) 
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+<0|C fc iCfcaCfc3Cfc4(t)T{^ f° f d 3 xdtd 3 x'dt"H 3 (x)H 3 (x')}\0} 

" J —oo J —oo 



to 



+{0\T{t / d s xdtH 4 (x)}( kl ( k2 ( k3 ( M (t)\0) 

J —oo 

+(0\tkitk2tk3tk4(t)T{-i I ° d 3 x dt"n 4 (x')}\0) + ... (13) 

J — oo 

where H3, T-L^ are the third and fourth-order Hamiltonian in the interaction picture. The 
latter two terms make up the contact interaction diagram, the rest is responsible for the 
scalar exchange. Let us also remind the reader that the gauge invariant observable 
( is, at first approximation, linearly related to the scalar tc via £ = —Hit. Also, 
already at this stage one can see that the result of the four point function is going to 
depend on six variables. All wavef unctions, once in Fourier space, depend only on the 
magnitude of their momenta. There are at most ten fields involved in the contractions, 
eight of which will always depend on the magnitude of the four external momenta 
(ki,k2,ks,ki). We are left with one last contraction between two fields depending on 
the magnitude of one vector which, by construction, is going to be the sum of two 
external momenta. It turns out that, employing the overall momentum conservation, 
two variables are sufficient to describe any of these linear combinations, we choose 
(A; 12 = \ki + k 2 1, ki4 = \k% + A^l), giving a total of six variables. As clear from above, 
the Mg-driven third-order interaction we are going to consider further depends on scalar 
products between the various momenta but, as one can easily verify, these can all be 
fully specified by using the six variables introduced above. All the variables we will 
employ are represented in the figure below. 



k 



2 



Figure B: the regular tetrahedron described by the four external momenta and k\2, fci4 



In order to get a tetrahedron as the one in Fig B one must enforce the following 
inequalities: 

cos(a — 0) > 005(7) ^ cos(a + (3) (14) 

with 

p±p _ p pip _ p p , p _ p 

cos(a) = fcl + \ cos^)= ^ + ^l fc3 ; cosfr) = 1 + * 2 , * 12 - (15) 

LK\K\At Z/C2/C14 LK\K2 
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From — 1 < cos(a, j3,j) < 1 one also obtains the usual triangles inequalities. Here we 
single out some of the inequalities which we are going to use in what follows: 

kj 2 + kj 4 < 4 ; ^l + kl-2fa < k 14 < ^l + k\ + 2k 4 . (16) 

In order to have a visual intuition and understanding of the result, once the calculation 
of the several contributions to the trispectrum is performed one needs to set up a number 
of configurations in which four out of the six variables are held fixed. Having more than 
one configuration also increases one's ability to distinguish the signatures of different 
interactions. Following [33], we adopt the set up described below: 

• Equilateral configuration: all the external momenta have the same magnitude 
k — k± — k 2 = &3 = k^ the two variables left are plotted as k\ 2 /k,k\ 4 /k. 
Note that when plotting in this configuration we will use the first inequality in 
Eq. (16). Incidentally, this is the only configuration for which exact calculations for 
the trispectrum in ghost inflation have been presented (see [30]) so far. Note also 
that for the equilateral as well as for the other configurations, one conveniently plots 
the result of the calculations in Eq. (13) for any specific interaction term multiplied 
by a factor of Y\t=i ■ ^ i s done also because this factor is generally common to 
all the contributions and so removing it sharpens the differences between the plots 
of each interaction term. 

• Folded configuration: here one has k\ 2 — > as well as k\ = k 2 and k 3 = k 4 . The 
second and third inequalities in Eq. (16) must be enforced in this case. The variables 
ku and k 4 are the ones plotted in this configuration. 

• Specialized planar limit configuration: in this case we have k± — k 3 — ku as well as: 



k 



12 



k\ + Y^ (^4 + ^(Akf-kl)(4kf-kl) )] . (17) 



2k\ 

The variables plotted are going to be k 2 /ki and k 4 /k 1 . 

Near double squeezed limit configuration: the tetrahedron is now a planar 
quadrangle and fc 3 = k 4 = k\ 2 . The region of interest is in particular the one 
for which fc 3 , k 4 , k\ 2 — > where the following relation holds: 

, /h 2 ( — h 2 _|_ £-2 I L,2\ _pp T p p _|_ p UT p p _ p h.2 _ L.4 _|_ u2l,2 
y r*j\ \ ">\ 2 i >"3 i ro 4 ) ™sl s2 ' 12 14 ' r°l2^4 ' "'14"'4 " / 14"'3 "'4 ' rv%rv 4 



k 2 = 



y/2k 4 



(18) 



with 

k 2 sl = 2^(k 1 k 4 + ki • k4)(fcifc 4 - k x • k 4 ) k 2 s2 = ^/(k 3 k 4 + k 3 • k 4 )(A; 3 A;4 - k 3 • k 4 ). 

(19) 

In this well the last two inequalities of Eq. (16) will be imposed on the 

variables k 14 /ki and k 4 jk\. Note that, only in this configuration, what one actually 
plots is the result of Eq. (13) times n£=i instead of nLi ^f- This is once again 
done in order to better appreciate the difference among the many interaction terms. 
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We now consider the result for the scalar exchange contribution focusing in particular 
on an interaction term (the one proportional to M 6 in Eq. 6 ) case which proved very 
interesting in plotting the shape of the bispectrum [28]. 

In all the calculations that follow we use a simplifying assumption which has been 
verified to hold for 3-point functions and is expected to hold for higher correlators as 
well [28]. Instead of using the generalized wavefunction which comprises the general 
single-field inflation solution, the Ghost inflation one, etc. as its simplified limits, we 
employ the usual Hankel function H 3 / 2 {kc s r) as a solution to the equation of motion 
for the quadratic action. The rationale for such a simplification is that, as one can 
readily verify, the main contribution to higher order correlators comes as usual from 
the horizon-crossing region and precisely in that region the behaviour of the general 
solution of [28] (see also [41]) resembles very closely the one of the simpler specific DBI 
wavefunction. We elaborate further on this fact in Appendix B where some examples 
and comparisons of explicit calculations are provided. 

Before moving to the detailed analysis of the shape-functions for several interactions 
terms in various configurations, let us comment briefly on the amplitudes generically 
associated with these interaction terms. As noted before, building on the freedom on the 
M n coefficients allowed by the theory and on the possibility of employing a small speed 
of sound, Cj C 1 (the same holds for the parameters which represent the generalization 
of Cg, i.e. cto and /3q), one can obtain large values for the amplitude associated to each one 
of the curvature-generated terms we are going to study. This has been quantitatively 
verified for all the terms of Eq. (6) in [28]. As an example, consider the Mg-driven 
fourthorder interaction term. To estimate the size of the amplitude associated with a 
given interaction term one considers its ratio with the quadratic terms of the theory at 
freezing [24]. Applying this prescription to our example one obtains: 

£j&2(w)'(v» w )' M|/J 6 7r 4 /c s 6 M 2 H 2 ( 2 



£ 2 M| if 2 7r 2 M 2 4 c~ e 



(20) 



where the linear relation ( = —Hn has been used; taking ( ~ 10 5 gives a rough 
estimate of the size of the non linear corrections. 



4-2. Scalar exchange diagram 

Here we are going to consider the interaction term M 6 Tr(dijir) 2 /a 4 . Note that, as 
opposed to the M 8 , M 9 -regulated terms in Eq. (6) which give a flat shape for the 
bispectrum much like the M 6 -driven interaction, this term is actually invariant under 
the symmetry S2 while for it to be (approximately) invariant under SI one needs 
to require its M$ coefficient to be much smaller than M4 which would in turn make 
its signal undetectable. We now write more explicitly the contribution of the Mq- 
driven third order interaction to the scalar exchange diagram. For all the details of the 
calculation, including contractions, we refer the reader to Appendix A. Consider here 
just one particular contraction of the fields, the sample contribution we are after looks 
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like the following: 

< n kl n k2 ir k3 ir k4 >Jg£= 

a// contractions 



3 • _ _ /i, i, n2 /„4 



n-,.,7r;../T/,..7ri.,(U) I dtia^TTfcia^iTTfcaCki ■ k 2 ) i /a 4 / (it 2 a 3 7rfc 12 7r^ 3 7r^ 4 (k 3 ■ k 4 ) 2 /a 4 



(21) 



-2K 



^* 



t->0 



d *i^fci fl 7r fci 7r fa(*i)( k i ■ k 2 



^^^a 71 ** 71 *^)^ • k 4 ) 



(22) 



where an overall momentum conservation delta and a factor of (27r) 3 have been omitted 
for simplicity. 

In Fig. 1, 2 below we plot what one obtains by summing over all contractions, accounting 
for the symmetry factors of the vertices and plotting the result. 

For the sake of comparison we often make reference to the shapes obtained in [33] and 
[30]. In the former work the so called local trispectrum is also plotted and compared 
with findings for general single-field inflation models; in the latter one shapes for 
the trispectrum of ghost inflation are presented in the equilateral configuration only. 
As mentioned, our starting Lagrangian comprises both these inflationary models; we 
decided to concentrate on plotting the novel curvature-generated terms that are invariant 
under SI, S2 or both. 





Figure 1. The equilateral configuration shape is presented on the left from the scalar 
exchange contribution of the Mg-driven interaction. The shape-function is different 
from the plots presented in [33]. On the other hand, once a necessary change of 
variables has been performed, it is qualitatively similar to the shape for the contact 
interaction diagram which arises from the ghost inflation (V7r) 4 interaction term in 
[30]. 

On the right we plotted our findings for the Mg-generated interaction in the folded 
configuration. It very much resembles the ones obtained in [33] for the scalar exchange 
diagrams from DBI-like terms, especially from the interaction 7r(V7r) 2 . In all the 
pictures above and below k\ has ben set equal to unity without loss of generality. 
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Figure 2. On the left the plot obtained for the Mq interaction in the planar 
configuration plot. Note here some of the interesting features: as &2,&4 — > the 
shape function goes to zero, much like it happens for DBI-generated interactions. As 
k2,ki — > 2 the shape function reaches values that are negative, albeit slightly so. 
Looking at the ki = k^ line we see that it is convex, rather than concave as found for 
other interaction types in [33] 

On the right the planar limit double squeezed configuration is plotted. In the region of 
interest, namely for k\2 — > 0, the shape function is non-zero, finite and negative; this 
again is different than what found in [33] for a variety of DBI-originated terms. 

Overall we see that performing the shape analysis for this S2-abiding term we are 
able to find some distinctive features with respect to the DBI-generated contributions. 
This is quite interesting also given the fact that the very same M6-modulated interac- 
tion term gives rise to a flat shape for the three point function [28] which contributes 
to enlarge the allowed classification of bispectra shape-functions for single-field models 
of inflation (with Bunch-Davies vacuum). 



4-3. Contact interaction diagram 

We now turn to the calculation of various terms that contribute to the contact interaction 
diagrams. The terms driven by M 3 , M4 are DBI-generated and have been calculated in 
a number of papers, notably [16, 33]. M2 is found in both DBI and Ghost inflationary 
theories [33, 30]. If we are to preserve both symmetries then we need to focus on 
Mil, M12, N1..N5 as one can easily check from Table 1. 

Given any interaction term at fourth order, its contribution to the contact interaction 
diagram can be written as: 



( Vk2Vki (t -> 0))£ = -2X m {0\T{i / d'x dt H 4 (x)7r* 1 7r* 2 7r* 3 7r* 4 (t -> 0)}|0) 

J —00 

(23) 

where X m stands for taking the imaginary part; a delta enforcing momentum conserva- 
tion and unimportant numerical factors have been omitted. 
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We have applied the above formula to a number of fourth-order interaction terms provid- 
ing some examples of notable SI and S2-abiding terms, SI or S2 invariant contribu- 
tions and, finally, terms that do not respect any of the symmetries above. We start with 
the S2-invariant interaction term ~ M|(V7r) 4 . This contribution is present in both DBI- 
like and Ghost inflationary models (e.g. in DBI one simply has l/c 2 s = 1 — 2M^ / MpH) . 
Although already written down in [30, 33], this term has not been plotted in all four 
configurations described in section 4 we employ here. This is because in DBI-like theo- 
ries it is expected to be subdomimant with respect to the M^-ir 4 term. 

• Ox = 1/2 M 4 (<9 i7 r) 4 /a 4 

For the shape function of the operator (V7r) 4 we see that a number of interesting issues 
arise. First, the plot in the equilateral configuration does not resemble any of those 
plotted in [33]. * Then, as we mentioned in Fig. 2, in the double squeezed configuration 
the k\2 — > limit gives a non-zero finite shape function. This is important because, up 
to the results in [33], this limit was thought as very useful to distinguish the leading 
contributions coming from interactions at third order in perturbations from the ones at 
fourth order in fluctuations. 



2.0 




Figure 3. The equilateral configuration shape for the 0\ operator is presented on 
the left. It is different from the results in [33] and very much resembles the plot we 
obtained for the scalar exchange calculation. Notice that upon performing a change 
of variables, it is basically identical to the (V7r) 4 interaction term plotted in [30] (We 
elaborate further on this point in Appendix B). 

On the right we plotted our findings for the 0\ interaction term in the folded 
configuration. As it will be for the other interactions, this configurations provides 
no particularly distinctive features that would allow to single out the constributions 
from the different interaction operators. 

In fact, all of the terms contributing to the scalar exchange diagram in [33] give a shape 
function which in the k\2 — > is finite. On the contrary, the leading contact interaction 

* See Appendix B for a detailed account of the plot of this term first done in [30] with different variables 
giving results which are essentially identical to ours despite a simplyfing assumption on our part. 
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diagram contributions analyzed in [33] do vanish in this limit. Note also that if one is 
to relax the assumption of a Bunch-Davies vacuum, the authors of [33] showed that this 
is not true anymore. 





Figure 4. The 0\ interaction planar configuration shape on the left is not exceedingly 
different from the ones presented in [33]: it vanishes for ^2,^4 — > 0, it is peaked for 
&2 = 2 = k±. On the &2 = &4 line the shapefunction is convex, rather than concave as 
for the contact interaction term plotted in [33] . 

On the right we plotted the 0\ interaction shape function in the planar limit double 
squeezed configuration. Here we immediately note an interesting feature: despite this 
being a contribution to the contact interaction diagram, in the = ki2 — > it gives 
a finite, non zero shape function. We comment more on this fact in the text. 

Next, we continue keeping our attention focused on terms which are SI and S2 in- 
variant. These include, in terms of their free coefficient, Mu, M12, N1..N5. Since they 
generate shape functions which are qualitatively very similar, we chose to plot just two 
representative terms in this list. 

• 2 = l/6 M^n 2 (dfn) 2 /a 4 



In the plots below we see that the Mn-driven term, generates shapes which are 
very similar to the ones plotted in [33] for the DBI-generated term ~ 7r 4 . 
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Figure 5. The equilateral configuration shape is presented on the left for the O2 
operator. Because of the way the space derivatives are written in Fourier space there's 
no fci2 , fci4 dependence and so one gets a plateau. On the right our findings for O2 in 
the folded configuration. 




Figure 6. The O2 planar configuration shape on the left: it vanishes for k-z, &4 — » 0, 
it is peaked for k 2 = 2 = k^. On the k 2 = k^ line the shapefunction is now concave, 
just like one would get for the tt 4 interaction. 

On the right we plotted the shape function in the planar limit double squeezed 
configuration associated to the O2 interaction term. 

We proceed with the other representative term: 

• 3 = 1/4! NsdfrdijTrdjkTrduTr /a 8 

The differences with respect to the Mn-driven interaction shapes are to be find 
in the first and the third configuration: in the first configuration they are due to the 
^-dependence of the interaction, on the third configuration A^3 gives a plot similar to 
the one tuned by the M 2 coefficient (see Fig. 4). 
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Figure 7. The equilateral configuration shape for O3 is presented on the left. On the 
right our findings for the O3 operator in the folded configuration. 





Figure 8. The planar configuration shape on the left for O3. On the right we plotted 
the shape function of the O3 interaction term in the planar limit double squeezed 
configuration. 

We now turn our attention onto terms which violate one of the symmetries, SI in this 
case. Indeed, we analyze the interaction ~ M|Mg 7r 2 (<9j7r) 2 /a 2 



4 = (M 2 4 M|)/(2M 2 4 + M 2 P eH 2 - 3M?H) x ii 2 (d^) 2 /c 
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Figure 9. The equilateral configuration shape for O4 is presented on the left. On 
the right our findings for the folded configuration associated to O4. Both are very 
similar to we obtained for the Mu in Fig. 5 and to what was found for the M4 -driven 
interaction in the literature [33] . 

For this interaction term we see the interesting feature presents itself in the fourth 
configuration where the k\2 — > limit gives a finite shape function. 




2.0 0.0 



Figure 10. The planar configuration shape for O4 on the left. On the right we plotted 
the 04-generated shape function in the planar limit double squeezed configuration which 
gives again a non-zero and finite shape function even for fourth-order interactions such 
as the one under scrutiny here. 

We now proceed to plot our findings for one more term, precisely the leading fourth- 
order interaction term among the ones driven by M 6 . It is clear that whenever this term 
gives a leading third order contribution (something one can achieve given the freedom 
on most M ra 's ), it violates Si. S2 however, is preserved by the leading terms associated 
to this coefficient at third and fourth order as is clear from Eq. (9) and Table 1. 

.0 5 = l/6M 6 (9 fc 7r) 2 (^7r) 2 /a 6 
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Figure 11. The equilateral configuration shape for the operator C 5 is presented on 
the left: this shape for the plot function has not been seen before in the equilateral 
configuration. 

On the right our findings for the O5 interaction term in the folded configuration. 



As one can see from Fig. 11, the plot in the equilateral configuration has no analogue 
in the shapes of [33, 30] for this configuration. It is somewhat reminiscent of the shape 
obtained for the 7r 4 of [33] but again, we stress it was obtained in a different configuration. 
The results plotted in Fig. 12 show once again that it is not safe in theories more general 
than DBI to attribuite to the planar limit double squeezed configuration the role to 
provide a distinctive signature in the ki 2 — > limit that would enable one to distinguish 
between third and fourth-order interaction contributions (see the discussion in [33]). 




Figure 12. The planar configuration shape for O5 on the left. On the right we plotted 
the 05-generated shape function in the planar limit double squeezed configuration 
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5. Conclusions 



In this work we aimed at employing the tools of effective field theory [24] to probe the 
space of all possible interactions for inflation driven by a single scalar degree of freedom 
up to fourth order in perturbations. Following a companion work on the bispectrum 
[28], and guided by the requirement of some additional symmetries on the action as 
an ordering principle [30, 31], in Section 4 we singled out the distinctive features in 
the trispectrum one obtains when considering curvature-generated terms from a very 
general fourth-order Hamiltonian. It is important to note that all of these intereactions 
allow, by construction, for a large trispectrum. Some of them present features which 
also emerge in DBI-inflation and Ghost inflation [33, 30]. 

We have also identified third and fourth-order interaction terms which have not been 
analyzed before and whose analysis reveals novel interesting effects. We were able 
to show that, unlike what happens in DBI-like models, the analysis of the double 
squeezed configuration cannot give a clear cut clue as to what kind of signal comes 
from leading third-order terms as opposed to fourth-order terms in perturbations. We 
found many interactions generating a shape in the equilateral configuration that mimics 
the behaviour of the ghost interaction term ( i.e. (V7r) 4 ) shape function first plotted 
in [30] which is quite different from the shapes of the DBI model (we also extended the 
ghost inflation plots of [30] to three other configurations). Triggered by an interesting 
flat shape found in [28], specifically the one controlled by the M 6 coefficient, the analysis 
of this term contributions to the scalar exchange and contact interaction diagram was 
calculated and plotted: a shape function which has not been found before emerged in 
the equilateral configuration for the contact interaction contribution. Finally, in Section 
4, we again verified up to fourth order a feature which already emerged at the level of the 
bispectrum in [28] : a shape function which in general single field inflation models is only 
obtained either by employing a linear combination of operators (as far as the Bispectrum 
is concerned) or relaxing the Bunch-Davies vacuum requirement for the theory, quite 
naturally arises in more general setups as the one employed here. Furthermore, it does 
so when considering several and independent interaction terms. 
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6. Appendix A 

A more detailed presentation of the scalar exchange calculation for the Mg-driven term 
is presented here. We start from the scalar exchange part of Eq. (13): 

(^ICfciCfc 2 Cfc 3 Cfc4WI^)s.e. = 

(0\f{i P d 3 xdtU 3 (x)}( kl ( k2 ( k3 ( k4 (t)T{-t P d 3 x'dt"U 3 (x)}\0) 

J — OO J —oo 



-OO 

1 1 j3„j +J 3, 



+ (0|T{- J J d 3 xdtd 6 x dt H 3 (x)H3(x)}( kl ( k2 ( k3 ( k4 (t)\ty 
HO\(kMk 3 (kAt)T{^ P P d 3 xdtd 3 x'dt"H 3 (x)H 3 (x')}\0}. 

" J — oo J —oo 

(24) 

Using Wick contraction on a generic operator 0, one has: 

T{0(ti)0(t 2 )0(t 3 )0(t 4 )} = ^{0(^1)0(^2)0(^3)0(^4) + all contractions}, (25) 

where N is the normal ordering operator. Note also that, being our 7r ~ ( operators 
squeezed between two vacua of the free theory, this reduces to considering only terms 
which are writeable as fully contracted contributions. For the anti-time order operator 
the same formula holds, only one needs to define contractions differently. We show 
below this difference: 



1 1 

0(x 1 ,t 1 ) 0(f 2 ,t 2 ) T = [0 + (fl,tl),0'(f 2 ,t 2 )]^(tl ~t 2 ) + [0 + O? 2 ,t 2 ),0-(fi,ti)Mt 2 -h) 

(26) 



1 1 

0O?i,ti) 0(f 2 ,t 2 )f = [0 + (:r 2 ,t 2 ),0"(:ri,ti)]6>(ti -t 2 ) + [0 + (:ri,ti),0~(:r 2 ,t 2 )]6»(t 2 -t x ) 

(27) 

where 

/d^k - f d 3 k 

— 0(M 2 )a k e^; 0"(f 2 ,t 2 ) = J j^4>*{K h) a£ e -* a . (28) 

Using the definitions above one gets several different contributions from Eq. (24). Note 
also that, using time and anti-time order definitions, the last two lines of Eq. (24) are 
just each other's conjugate and can therefore be grouped together. We now procede 
to write an explicit expression for the four point function generated by the M 6 -driven 
contribution to the scalar exchange diagram: 

(7rfcl7r fc2 7T fc 37T fc 4)^ 6 oc (29) 

„3 pO „3 

(0) / 

+ 5 permutations) + 



Mf 



( f° ° 3 2 f° ° 3 * 2 

4 • I ^1^2^3^4(0) / (iti— 7r fc i 2 7r fe i7r fe2 (ki • k 2 ) 2 / dt 2 —n* kl2 'K* k3 7r* k4: (k 3 -k i ) 2 
V J -00 a J— 00 a 
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( f fl3 2 f a 3 

12 ■ 2 ■ n* kl 7r* k2 ir k3 ir k4 (0) / dti— 7rfci7r fe i 2 7r fc2 (-ki 2 • k 2 ) 2 / ^2—^12^3^4(^3 • k 4 

+ 11 permutations) 

( f° a 3 . 2 f° a 3 . „ 2 

2 • 2 • 7^7^27^37^4(0) / rfti— 7r fc i 2 7r fe i7r fc2 (ki • k 2 ) 2 / ^2^7r^ 3 7r^ 12 7r^ 4 (ki 2 • k 4 ) 2 
V " J-00 a J-00 a 4 

+ 11 permutations) 

f° a 3 . 2 Z" a 3 . . . » 2 

4 ■ 7r^ 1 7r^27r fe 37r fe4 (0) / citi— 7r fc i7r fe i 2 7r fe2 (-ki 2 ■ k 2 ) 2 / ^ 2 ^-7r^ 3 7rfc 12 7r^ 4 (ki 2 • k 4 ) 2 

J-oo a J-oo a 



+ 23 permutations)} + 

pO „3 rti ^3 



2M| 



3 

+ 5 permutations) + 



4 ' ( 7r fci 7r fc2 7r fc3 7r fc4(°) / d *i~4 ^12^1^2 (ki • k 2 ) 2 / dt 2 — 7r A .i 2 7r A , 3 7r A , 4 (k 3 
V a 4 a 4 



/ f° a 3 /•*! a 3 

2 • 2 • 7r^ 2 7r^ 4 (0) / rftx— 7V fcl 7r* 12 7r fe2 (-k 12 -k 2 ) 2 / d* 2 — 7r fcl2 7r fe3 7r fc4 (k 3 • k 4 )^ 

V J-oo a J-oo a 

+ 11 permutations) 

2 • 2 • Tfci"'fc2 7r fc3 7r fc4( ) / dt i— ^12^1^2 (ki • k 2 ) 2 / cit 2 — 7r fc3 7r fcl2 7r fc4 (ki 2 • k 4 ) 2 

a J-00 a 



—00 j —00 

+ 11 permutations) 

„3 rii „3 



/" a 3 a 3 
4 • 7r fci 7r fc2 7r fc3 7r fc4(°) / 7 Tfci7rfc 12 7r fe2 (-k 12 • k 2 ) 2 / dt 2 — 7r fc3 7r fcl2 7r fe4 (k 12 • k 4 ) 2 

J-00 a J —00 a 

+ 23 permutations)} . 

One then performs these calculations and plots the results to obtain Fig. 1,2. The 
situation for the contact interaction diagram contributions is considerably simpler as 
there is just one time intergral to be performed and two less fields to be taken into 
account. 



7. Appendix B 

We want here to show with an example what seems to be a general feature concerning 
the use of (reasonably) approximated wavefunctions in the calculation of higher order 
correlators. In [28] we found that in performing an exact calculation for correlators 
in a very general theory such as the one we employed in this paper, whenever a 
given interaction term was producing a shape function for the trispectrum which one 
could qualitatively classify as, say, equilateral, so was the calculation performed with a 
simplified wavefunction. This is due to two independent reasons. First, we start from 
the realization that, precisely in the horizon-crossing region, which is where one expects 
the main contribution to any n-point to come from, the exact general wavefunction [28] 
and the usual one, H^l(c s kr), which in these theories is an approximated solution, 
behave very similarly. Secondly, in [28] we concluded that most of the distinctive effects 
of the bispectrum where due not to the particular k-modes dependence of the result 
of the integrals like the one in Eq. (22), but on the fraction of that k-dependence that 
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could be taken outside the integral, so on the part of the k-dependence not directly 
attached to the time behaviour of the wavefunction and which is common to the exact 
and approximated wavefunction. 

We now compare the trispectrum shapefunction of a ghost inflation interaction term, 
(V7r) 4 , performed with the exact ghost solution in [30] with the results we obtain 
employing the approximated DBI wavefunction, just what we used in obtaining all 
the shape functions presented here. 
The interaction reads: 



Its trispectrum shapefunction obtained through the approximated methods has been 
ploted in Fig. 3. In order to compare it with the exact calculation of [30] we need to 
change variables and turn to: 



C 2 = h-h; C 3 = h-h; C 4 = -1 - C 2 - C 3 ; k 12 = ^2(1 + C 2 ) k u = y/2-C 2 - C 3 ) 

(32) 



We now show the plots obtained by performing this change of variable on our 
approximated result alongside the plot obtained with the exact ghost wavefunction 
taken directly from [30]. 




(31) 



i=l j=l 



1.0 




1.0 



Figure 13. On the left the approximated result. The two shapefunctions are 
qualitatively identical barring an unimportant numerical coefficient due to a different 
normalization. 
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